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Abstract 



We apply the technique of self-similar exponential approximants based on 
successive truncations of continued exponentials to reconstruct functional laws 
of the quasi-exponential class from the knowledge of only a few terms of their 
power series. Comparison with the standard Pade approximants shows that, 
in general, the self-similar exponential approximants provide significantly bet- 
ter reconstructions. 
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I. INTRODUCTION 



Exponential laws are ubiquitous in nature. The overwhelming majority of relaxation 
phenomena occur through exponential laws. The dominant role played by exponential laws 
results from a combination of mechanisms. First, exponential relaxations result from viscous 
dissipation proportional to the first-order time derivative of the dynamical quantity, which 
is usually the dominant term breaking the time-reversal symmetry (see however |I| for 
a case where the third-order derivative becomes important). Second, exponential decay or 
growth reflects the first-order expansion in the rate of change as a function of the observable. 
Third, exponential laws are often associated with the Poisson process, which has the unique 
property of being memoryless. This leads to a remarkable mathematical property: the 
Poisson law is invariant with respect to any conditioning on the past or future and can thus 
be seen as the unique fixed-point of arbitrary transformations of time distributions involving 
conditioning 0. 

These universal properties provide useful benchmarks against which deviations can be 
gauged. These departures usually betray specific mechanisms at work in real systems, that 
are otherwise hidden by the just-mentioned universal mechanisms. This is why the study 
of such laws which are close to exponentials is so important. It may happen, however, that 
we are not able to observe in some experiment or to study theoretically such a law, say a 
relaxation process, from its beginning till its very end. In theoretical works, this repeat- 
edly happens when considering complicated problems that do not allow for exact solutions. 
Then, invoking some kind of perturbation theory, one may calculate a few successive ap- 
proximations. The standard situation is when one employs a short-time expansion resulting 
in approximations presented as polynomials over time. The main problem in such cases is 
how to reconstruct the overall process from the knowledge of only its short-time behavior 
represented by a series in powers of time. An analogous situation may also arise in exper- 
iments, when one cannot, because of some technical difficulties, continuously measure the 
whole temporal process, but one is able to gain information only from a limited number of 
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measurements performed at discrete times. In that case again, one often presents the results 
in the simplest form of a polynomial describing the given set of discrete points. 

It may also happen that the whole relaxation process is impossible to observe experi- 
mentally because the lifetime of a system is shorter than the relaxation time of the process 
being studied. For instance, intensive investigations of various physical processes are now 
being accomplished for trapped atoms (see reviews [^). Since the life time of atoms in 
a trap is finite, not all processes can be observed in full, in particular phenomena associated 
with trapped spinor Bose condensates @J^. 

One more example of the same problem arises when the characteristic relaxation times 
are too long to allow for convenient experiments covering the full range of the relaxation. 
Let us mention the spin-lattice relaxation in polarized nuclear magnets [P,^], in which the 
relaxation time may attain several days at low temperatures. Another illustration is the 



magnetic relaxation in molecular magnets |jTO[- IT2|) where the relaxation time ranges up 



to several months below the blocking temperature. There are cases when the characteristic 
times become comparable to or larger than human time scales, as for instance in visco-elastic 
relaxation in the earth crust after earthquakes |]T^,|T3 • Here, seismologists may only observe 



the early part of the relaxation process which may last for decades to centuries. In general, 
if one extracts information only from the beginning of a process, the result can again be 
presented as a power series. And the question that arises is how to find the general relaxation 
law from the knowledge of only its short-term presentation. 

Here, we suggest a technique to address this question of how to derive the whole sought 
function based on a perturbative expansion for short-time dynamics. To this end, we exploit 
the remarkable property of self-similar exponential approximants introduced in Ref. ||15|| , 



which is able to reconstruct exactly the exponential function exp(— t) from the knowledge of 
only the two first terms of its Taylor expansion. In other words, the technique of self-similar 
exponential approximants represents an ideal filter for the exponential function. It is thus 
natural to expect that this technique can be successfully applied for the reconstruction of 
continuous functions deviating from a pure exponential while still keeping the same expo- 



nential asymptotics at infinity, from tlie sole knowledge of their Taylor expansion for short 
times. We note that the self-similar exponential approximants generalize and give a sys- 
tematic justification of the method of summation of power series by continued exponentials 



proposed by Bender and Vinson |T^ as an alternative to Pade techniques, which themselves 
the result of truncated continued fraction representations of power series. 

In the next section, we briefly review the main properties of the technique of self-similar 
exponential approximants that will be useful for our purpose. In section 3, we then discuss 
several examples illustrating its power to reconstruct quasi-exponential laws. In section 4, 
we apply the technique to an ordinary nonlinear differential equation motivated by a solid- 
friction problem. In each case, we compare our results with those from the more standard 
Pade approximants. In general, we find that the numerical errors of the reconstructions are 
significantly smaller for self-similar exponential approximants than for the standard Pade 
approximants. 

II. SELF-SIMILAR EXPONENTIAL APPROXIMANTS 



The complete mathematical foundation of the method can be found in Refs. [|T^, |T7 



| 2l| . Here, we give a brief sketch which emphasizes the variant that is mostly appropriate 



for filtering exponential-type laws. Assume that we are interested in a function (j){t) of a 
real variable t. Let perturbation theory (or some fitting procedure) give for this function 
the perturbative approximations (f)n{t), with n = 0, 1,2, ... enumerating the approximation 
order. Consider the case when the perturbative procedure results in a polynomial 

n 

0„(t)~^afct'= (t^O), n = 0,l,... (1) 

k=0 

Let us stress that the expansion (1) has, as a rule, no direct meaning if continued straightfor- 
wardly to the region of finite arbitrary t. In theoretical physics, the problem of reconstructing 
the value of a function at some distant moment of time from the knowledge of its asymptotic 
expansion as t ^ 0, is called the renormalization or resummation problem. An analytical 
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tool for the solution of this problem, called algebraic self-similar renormalization, has been 
recently developed |T5[, ||T^- The polynomial representation (1) gives for the sought 
function the following n polynomial approximations i = 0, 2, n, 

Po(i) = flo, pi{t) = aQ + ait, p2(t) = pi(t) + a2t^ , . . . , p„(t) = p„__i(t) + a„r. (2) 

The algebraic self-similar renormalization starts by applying to the approximations (|^) an 
algebraic transformation, thus defining a new sequence, Piit, s) = t^pi(t), i = 0,2, n, with 
s > 0. This transformation raises the powers of the series (1), (2), and allows us to take 
effectively into consideration a longer timespan of the system history. We shall use below the 
strongest form of such transformation occurring by taking formally the limit s — *• oo. It can 
be shown that this results in an exponential representation of the sought function. The next 
step of the self-similar approximation theory [^]- consists in considering the sequence of 
transformed approximations, Pi{t, s), as a dynamical system in discrete time i = 0, 1, ...,n, 
that we call "time-order" to distinguish it from the time variable t. In order to define 
the system evolution in time-order, it is convenient to introduce a new variable (p and to 
define the so-called expansion function t{ip, s) from the equation Po(^) s) = a^t^ = if, which 
gives t{(p, s) = {Lp/ao)^^^ . This makes it possible to construct the cascade of approximations 
yi{{p,s) = Pi{ t{{p,s),s). Embedding this cascade into a continuous approximation flow, 
one can write the evolution equation in terms of the discrete time-order variable in the 
form of the functional self-similarity relation, yi^p{ip,s) = yi{yp{ip, s), s), which is also the 
necessary condition for the convergence of Pj. Already at this stage, we can try to check the 
effectiveness of the algebraic transformation by calculating the so-called local multipliers, 

(3) 

f=Po(t,s) 

as s oo. When all \mi{t, oo)\ < 1, the convergence of the sequence Pi is guaranteed. To 
implement the calculations concretely, one can use the integral form of the self-similarity 
relation, 

i-p: dip 
Jp,_iVi{p,s) 
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mi{t, s) = 



dyiip-s) 
dp 



where the cascade velocity is Vi{ip,s) = yi{ip,s) — yi^i{ip,s) and r is the minimal number 
of steps of the approximation procedure needed to reach the fixed point P*{t,s) of the 
approximation cascade. It is possible to find P*{t,s) explicitly and to perform an inverse 
algebraic transform after which the limit s ^ cxd is to be taken. The first step of the self- 
similar renormalization is completed. This procedure is then repeated as many times as 
necessary to renormalize all polynomials which appear at the preceding steps. Completing 
this program, we come to the following sequence of self-similar exponential approximations 

Pj{t,TiT2,...,Tj) = floexp ^^ritexp ^^r2t...exp ^^^^i^^^ •••^ ' J = 2,3...,n. (4) 

The quantities ti, T2, . . . play the role of control functions that are to be defined from 
optimization conditions. For the purpose of reconstructing an exponential-type law, the 
optimal choice of controls is obtained by expanding p*{t,Ti,T2, ...,Tj) in the vicinity of 
t = and demanding that this expansion coincides with (pjit). In the theory of Fade 
approximants, this condition is often refered to as the "accuracy-through-order relationship." 
This method of determination of the control coefficients Ti makes the self-similar exponential 
approximants analogous to the Euler superexponentials |lT6| , p6| , |27[| . With this determination 



of the control parameters, we come finally to the self-similar approximants 0* of the sought 
function, 

(l)*j{t, T1T2, ...,Tj)= p*{t, Ti, Ta, . . . , Tj) . 

For example, 

02(t, n, T2) = ao exp ( —Tit exp ( — rgt 

(j)l{t,Ti,T2,T3) = aoexp (— Tit exp (—Tat exp ( — rgt 

Veto Vfli \a2 

( a\ ( a2 / CI3 / 04 

(/)4()f:, Ti, r2, Ts, r4) = aoexp — ritexp —rat exp — rat exp — r4t 

Voo Vai \a2 V0.3 . 

In order to check whether the sequence of 0*(t, ri, ra, Tj) converges, we study their map- 
ping multipliers, M*{t, ti,T2, . . .Tj) defined as 

A 6^*{t,n,...,Tj) 1 d 
M, . . . ,r,) ^ -^-S^) = ^ ^ • • • '-^•) • 
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This definition of tfie multipliers allows us to compare the convergence of the expansion and 
of the renormalized expressions, making clear what can be expected a priori. 

With the control parameters defined as prescribed above for each j from the the accuracy- 
through-order relationship, we can obtain j self-similar exponential approximants for the 
sought function (where all r are now known functions of the parameters Oj): 

0*1 (t) = <^*(t, n, 1, 1), (j)*,{t) = (j)*{t, n, r2, 1, 1), = (j)*{t, n, r^, r,), 

which differ according to the number of control parameters being employed. This provides 
a matrix of self-similar approximants, indexed by the order j and by the number of control 
parameters. To this matrix of approximants is associated the matrix of multipliers. E.g., 
for j = 4, we have 

<Pli{t) = <Pl{t, ri, 1, 1, 1) , <l)l,{t) = <pi{t, n, T2, 1, 1) , 

043 W = (l^li^^ ^1' ^2, T-3, 1) , (t) = (f)l{t, Ti, T2, Tg, n). 

We propose to qualify the convergence towards the sought function by examining both the 
convergence of the different sequences of multipliers M* and of the approximants 0*, choosing 
among them the pair which exhibits the best convergence rate. The resulting limiting point 
from the table of approximants should be taken for the value of the sought function. If there 



are more than one limiting point, one can take their weighted average, following Ref. |21 



III. ILLUSTRATION OF TECHNIQUE IN ACTION 

We now present explicit examples demonstrating how the technique works. 

A. Ideal Filter 



Consider the Taylor expansion of the exponential function, (j){t) = exp(— t), up to an 
arbitrary order in t, 
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where 



,(t)^5:a,t\ (6) 

fc=0 



ttQ — i, ai — — i, 02 — -, 03 — — -, 04 — — , 05 — — 



The second-order self-similar approximant is 

f CLl f 0,2 \\ la? 
4>22V^) = '^0 exp t Ti exp ^ 'i^2 , Tl = 1, T2 = 1 - - = 0. 



^ (2q veil / / ^ (X2CiQ 

The fact that r2 = leads to 022('^) ~ exp(— t). It also makes any arbitrary order approx- 
imant <t>*jj{t) identical to exp(— t). Note, that all other approximants, except 02i, ^31, ^41 
and 05]^, become identical to exp(— t) and all higher order control parameters are identical 
zeros, like T2. 

Consider now the conventional Pade approximants, P|(t) and P4(t) (Pi(t) and P2(^) have 
inferior quality) and compare them to the exact function, see Fig. 1. Although P|(t) remains 
positive, it behaves nonmonotically. On the other hand, P4(t) becomes negative for large 
t. These are typical problems encountered while attempting to reconstruct functions with 
exponential asjTiiptotic behavior by means of Pade approximants. The relative percentage 
errors for those two Pade approximants are shown in Fig. 2. One should not be mislead 
by the seemingly superior performance of P4(t), which is qualitatively wrong in predicting 
negative values already for moderate times. Because of this, we will present below only the 
relative percentage error for positively defined Pade approximants. 

Let us study the impact of noise on the coefficients of the power law expansion @. Since 
022 (^) recovers the exact solution in absence of noise, we study its perturbed value denoted 
022(^5 Vi ^) brought by the existence of the noises 77 and 9 on the two first coefficients of the 
expansion @ defined by the replacement of ai by ai(l + 77) and/or of 02 by 02(1 + ^^). For 
sufficiently small amplitudes of the fluctuations, one can expand 022(^5 Vi ^) follows (only 
the most simple cases are considered): 

cj>l^{t, V, 0) - 0;2(t) ^ hit) r, + hit) + 
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and 

<p*,{t, 0, 6) - <p;^{t) ~ ci(t) e + C2{t) e^ + ..., 

where all coefficients can be easily calculated analytically. We observe that: 

1. the absolute error in the lowest order tends to peak at intermediate times, see Fig 3, 
while the relative error monotonously increases with time; 

2. the absolute error caused by fluctuations in 02 is smaller than the error caused by 
fluctuations of equivalent magnitude in oi, see Fig 3; 

3. higher order terms in the expansion tend to shift the peak in the absolute error 
towards larger times, see Fig. 4; 

4. Interference of the error peaks at different orders leads to broaden the domain in time 
over which the errors are significant, see Fig. 4. 

It is possible to perform an averaging over the random variables and evaluate the cor- 
responding averages as a formal series to develop a specific resummation procedure that 
may provide a solution for the average function, even for arbitrary strong noises. Wc do not 
persue this here and return to the study of how self-similar exponential approximants, which 
are ideal filters for the exponential function, extract information on other functions which 
are perturbations of exponentials, from a few starting terms of their Taylor expansions. 

B. Weak perturbation 

Consider the function, = sin(i)^ exp(— i^), which is nonmonotonic. Its expansion is 

Mt)^t'j:o^kt'' , (7) 

fc=0 



where 



4 79 8 13921 29341 

ao — 1, Ql — — - ) ^2 — — , Ct3 — — — , 04 — o ^nn ' 0.5 — — 



3 ' 90 ' 21 ' 113400 ' 935550 

The coefficients alternate in sign, like in the ideal exponential example. The values of the 
control parameters, 

T2 = 0.013, Ts = -0.137, T4 = -0.131, T5 = -0.152, 



are all negative, with T2 being rather small, which guarantees a good quality of the approxi- 
mation already in the lowest nontrivial order. Higher-order approximants provide corrections 
to the good starting approximation 022(^)- 

Both multipliers and approximants converge extremely well. It is difficult to distinguish 
in Fig. 5 the approximant (pl^^it) from the exact function. The relative error at t = 2 is 
0.045% for this approximant, while (f)5{t = 2) defined in (0) has a very large error —4.058 x 
10^ %. Note that the nondiagonal forms and approximate the sought function rather 
badly compared to the diagonal approximants. The logarithms of the relative percentage 
errors for the Fade approximant P^t) and for the super-exponential approximant (j)l^{t), 
corresponding to the same number of terms used in the construction of both approximants, 
are shown in Fig. 6. 



C. "Badly damaged" expansion 

Let us consider a function with an expansion whose coefficients do not alternate in sign, 
such as (j){t) = ln(l + t) exp(— t^), which has the expansion 

05(t)~t^afet\ (8) 



fc=0 



with the coefficients 



_1 _2 _1 _11 _1 

z 6 A 60 

This is again a nonmonotonic function. The control parameters, 

T2 = 1.187, rg = 1.69, u = -0.296, rg = 0.661, 

have different signs, with T2 and being rather large, which leads to a poor quality of 
approximation to the exact function for 022 (^) ^^"^ ^'ssi't)- Rather than finding a good ap- 
proximation to the exact solution already at the ffist step of the construction, the successive 
self-similar exponential approximants bracket the exact function within a rather broad range, 
as seen from the change of the sign between and r4. The analysis of the multipliers shown 
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in Fig. 7 shows that Ml^{t) and M^^{t) are the closest, which suggests a better convergence 
and that (pl^it) and 044(t) are the two best approximants. Indeed, they are much closer 
to the exact function than the lower-order approximants as shown in Fig. 8. Also, since 
|M55(t)| < |M44(t)|, one can anticipate correctly that 4>lr^{t) is located closer to the exact 
solution than 044(t). 

Consider all Fade approximants {Pl(t), F4(t), Pi(t), and F|(t)) which can be built from 
the fifth-order polynomial One can see in Fig. 9 that their quality is by far inferior to 
that obtained with ^55. Fig. 10 shows the logarithm of the relative percentage errors for 
the self-similar exponential approximant (pl^it) and for the Fade approximant Pf{t), which 
is the only one remaining positive. 

D. Strong perturbation 

Consider the Taylor expansion ip^it) ~ J2k=oO'k of the exponential-type function, 
= i^exp(— t), up to the fifth order in t, with the coefficients 

_ 5 _ 8 _ 65 _ 163 

The coefficients alternate in sign but their amplitudes behave nonmonotically. The approx- 
imants (j)l^{t), (j)l^{t), <j)Ut), <j)Ut), and multipliers M*^{t), M*,{t), Ml^{t), M*^{t) for the 
sought function can be readily written down. The values of the control parameters, 

T2 = 0.2, Tg = 0.508, r4 = 0.377, rg = 0.378, 

are all positive and smaller than one. 

Fig. 11 shows the dependence of the multipliers as a function of time. One can observe 
that M44(t) and M^^(t) are the closest pair, suggesting a good convergence of the cascade of 
approximants. The corresponding approximants, 044(t) and (pl^it) bracket the exact function 
from above and below respectively, see Fig. 12. One can observe that the approximant (p^^it) 
has already a rather good quality of approximation of the exact function. Relative percentage 
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errors for the Pade approximant P|(t) and for the self-similar exponential approximant 4>l^{t) 
are shown in Fig. 13. In a narrow region, Pl{t) outperforms (pl^^t) but fails short for larger 
times. However, due to the fact that the control parameters T2, t^, and remain all 
positive and of similar magnitude, the quality of 055(t) becomes eventually worse than that 
of the Fade approximation at very large times. 

To investigate this behavior some more, let us construct the following function, = 
exp(— t (1 + 1)^^^^), which decays as exp(— ^/t), as time goes to infinity. Its Taylor expan- 
sion, 05 (it) ~ Y%=QO'k possesses the coefficients 

25 53 2261 

On = 1, Oi = —1, ao = 1, 03 = , 04 = — , as = . 

u , i ,2,6 24' 48' 1920 

The control parameters, 

Ta = 0.5, T3 = 0.48, T4 = 0.393, rg = 0.367, 

are all positive, smaller than one but decrease rather slowly. Fig. 14 shows the dependence 
of the multipliers as a function of time. One can observe their poor convergence. As a 
consequence, the self-similar exponential approximants (/>44(t) and (t) bracket the exact 
function rather poorly, see Fig. 15, with an accuracy inferior to that of the Fade approximant 
P\{t), with the exception of very large times, see Fig. 16. 

These two examples illustrate the property that, for coefficients Oj's rapidly growing 
in absolute values, the considered self-similar exponential approximants, with controls de- 
scribed by the accuracy-trough-order relationship, become unreliable. Such cases of coeffi- 
cients growing as fast as a factorial of their order constitute an important class of behavior, 
since it appears in expansions that are typical to many nonlinear field theories. Consider for 
instance the particularly illustrative example of the Stieltjes function, (f)(t) = /q°° ^^^^~_^^ du, 
which exemplifies such a behavior. Its coefficients of the Euler series 

ao = 1, ai = -1!, a2 = 2!, = - 3!, 04 = 4!, 05 = -5!, 

diverging as a factorial, lead to the control parameters 

Ts = 0.75, ra = 0.713, n = 0.697, = 0.685 . 
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As we can anticipate from our previous observations, both the multiphers and approximants 
sequences are not convergent simultaneously, as seen in Figs. 17 and 18. In this case, the 
Fade approximant Plit) easily outperforms (pl^it), see Fig. 19. 

We stress that this failure of the Euler-type self-similar exponential approximants, as 
compared to the Fade approximants, does not imply the failure of the algebraic self-similar 
renormalization as a whole. Rather, it demonstrates the limitations of the particular way 
of defining controls and also of the infinite-power condition (s oo) used in the derivation 
of these approximants. In applications, the usage of the exponential approximants is often 
analogous to mean-field-type approximations. Because of this, we may call, for brevity, the 
limit s — s> oo as the mean-field condition. It is important to realize that this mean-field 
theory sends clear warnings about its own anticipated failure, reflected in the divergent 
behavior of the approximants and of their multipliers. In a subsequent paper, we show that 
lifting of the infinite-power restriction improves dramatically the accuracy of the self-similar 
renormalization. 



IV. EXPONENTIAL APPROXIMANTS DERIVED FROM A DIFFERENTIAL 

EQUATION 

We end our exploration of illustrative examples by the analysis of an expansion derived 
from the time evolution equation of the state variable usually called 6 in the friction literature 
of a block subjected to constant shear over normal stress, given by the Ruina-Dieterich solid 



friction law (see p9[ page 283 section 13.6.2 and ||30|j3l|] ). Fosing x = O/Oq where Oq is a 
parameter of the constitutive friction law, the following equation 

^ = - a x'-"", x{0) =xo>0 , (9) 

describes the evolution of the state variable of the friction law. In the following, we shall 
consider the case when the condition 9q^ — axl~"^ < holds which ensures that the solid 
friction state variable tends to decrease from an initial large value. This law (^) results from 
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a velocity-dependent solid friction coefficient, the block velocity V being related to the state 
variable x by the relation V oc l/x™. a is a parameter of the constitutive velocity-dependent 
solid friction law. The exponent m is also a parameter dependent upon the physical nature 
of the solid contacts. Equation (^) with initial condition — asg^™ < thus describes the 
acceleration of a block in contact with a solid substrate pulled with a constant force which 
has been suddenly applied at t = 0. It has been applied to describe one of the possible 



regimes of a mountainous slope which can become transiently unstable |^ 
In dimensionless variables, X = x/x, T = t/t, where 



1 1 

x = {aeo)~, t = \l-mf^a—eS'-\ (10) 



Eq. (PD reads 



^ (|1 - m|) = 1 - X^-, X(0) = ^0 = |, 1 - X^'" < . (11) 

For m < 1, it is a nonlinear relaxation equation, with x being an equilibrium value reached 
asymptotically after an exponential decay characterized by the typical time t. For m < 1, 
in a long-time limit, we obtain an asymptotic solution in the form, 

x(t) ~x + Aiexp (-t/t) + ..., t ^ oo , (12) 

where the value of Ai remains unknown. A naive approach consists then in writing down a 
naive approximate solution to Eq. (^, xi^{t), based on such an asymptotic single relaxation 
time expression, 

XN(t) =x + {xq — x) exp {—t/t) . (13) 

In this ansatz, the amplitude Ai is now determined. 

The exact expansion Xk{t) up to order t^ for short times and arbitrary m can be obtained 
from Eq. (y). Here, we limit ourselves to expansion up to the fifth order, 

5 

X5(t) ^ 5]aX, t^O, (14) 

n=0 
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where in dimensional units, 



Oo = Xq, ai = — a Xg 02 = -a {m — 1) oiXq 



as = -(^('Ti — 1) —m al Xq"^ ^ + 2a2 Xq ' 
04 = — a(m — 1) (1 + m) m a'^ Xq™'~'^ — 6m 02 aiXg "^^""^ + 603 Xq' 



as = ^a{m - 1) [-m(l + m) (2 + m)alxQ ™ ^ + 12m(l + m)xQ ^aia2- 



5! 



-24mXo'""^a3ai - 12mXo'""^a2 + 24x0"" 04 



Using this expansion, we apply our technique to construct a self-similar exponential 
approximant using this fifth-order expansion For m = 0.85, Xq = 50, Fig. 20 illustrates 
how well the "exact" numerical solution is approximated by superexponential approximants 
in dimensionless units. The values of the control parameters, 

T2 = -2.02, Ts = 1.504, n = 0.848, n = 0.649, 

have different signs, with T2 and T3 being rather large, which leads to a limited quality of the 
approximations at different orders. We stress that these values for the model parameters 
correspond to a highly nonlinear and strongly out-of-equilibrium case, where the naive expo- 
nential expression xi\i{t), corresponding to a close-to-equilibrium linearized model is way-off 
the mark both in qualitative and quantitative sense, see Fig. 20. In order to guarantee that 
the solution saturates at the correct constant value at infinite time, we add and subtract x 
from the initial series, which produces a single modification in that only ag is replaced by 
ao — X in the generic expressions for the approximants. We thus plot 

x{t)=x + xUt), z = 2,3,...5, (15) 

divided by x as a function of the dimensionless time T. Note that such shifts are completely 
compatible with the framework of algebraic self-similar renormalization, see e.g. Ref. |^ . 



For the fifth-order approximant studied here, the relative error is largest around T = 12.57 
and amounts to 1%. 
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In order to compare the accuracy of the self-similar approximants with the Pade approx- 



imants 28], we consider the Pade approximant, 



3 ao + Alt + A2t^ + Bsxt^ 

which goes to x at infinity, with all unknown coefficients determined from the short time 
expansion (0). Fig. 21 shows that completely fails at intermediate times and even 

goes through singularity before returning to asymptotic value x at large times. 

In summary, in this case, the self-similar exponential approximant has been able to cap- 
ture a highly non-trivial departure from a pure exponential relaxation which is dominating 
the relaxation process over a large time span, conditioned on the fact that the relaxation 
becomes asymptotically an exponential at long times. This shows again the power of this 
resummation method to capture significant deviations from exponentials in functions that 
belong to the exponential class in an asymptotic sense. 



V. CONCLUDING REMARKS 

In conclusion, we have demonstrated how the technique of self-similar exponential ap- 
proximants makes it possible to reconstruct exponential-type functions, when only a few 
terms of their expansions are known. 

This technique can also be applied to functions with asymptotic behavior different from 
an exponential. In this case, by carefully examining the convergence of the multipliers and 
of the approximants, it is possible to construct an accurate approximation for the sought 
function, while staying within the limits of applicability of the mean-field regime. 

Further increase of accuracy will come from lifting the mean-field condition used in 
deriving the self-similar exponential approximants. The result of this approach will be 
presented elsewhere. 

The Pade approximants remain a valuable technique, but it has no much value in theories 
of relaxation and should be replaced by other techniques, possibly by the superexponentials 
presented here and their non mean-field extensions that we shall report in a future work. 
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Figure Captions 

Figure 1. The conventional Fade approximants, -P|(t) (dash-dot) and -P4 (t) (dash) com- 
pared with exp{—t) (solid). Approximant 4>li{t) (dot) is shown as well. 

Figure 2. The relative percentage errors for the P^{t) (solid) and P^t) (dash) Fade 
approximants, are shown. One should not be mislead by the seemingly superior performance 
of Pi(t), which is qualitatively wrong in predicting negative values already for moderate 
times. 

Figure 3. Illustration of the impact of "noise" rj (or 9) in the coefficients of the power 
law expansion, onto the accuracy of the exponential approximant 022; characterized by 
the absolute error, 022(^>^5O) ~ ^22(^) — ^i(^) V + ^^it) rf + ... and 022('^'O>^) ~ 022(^) — 
Ci{t) 9 -\- C2{t) 9"^^-. The absolute errors in the lowest order, 1 61 (i)| (solid) and Ci{t) (dash), 
are shown as functions of time ( 77 = ^ = 1). 

Figure 4. Demonstration of the influence of the noise introduced into the higher order 
terms in the expansion onto the absolute error dependence on time. The absolute errors in 
the lowest order, |&i(f)| 77 (solid), the next order contribution to the error, &2(i) ff (dashed), 
and their sum (dotted line) are shown for 77 = 0.1. Same notations as in Fig. 3. 

Figure 5. It is impossible to distinguish the approximant 055(t) (dash-dot) from the 
exact function, (j){t) = sin(t)^ exp(— t^), intended to be shown with dots. Approximants 
(f)li{t) (solid line), 0|i(i) (dash line) and Taylor expansions 055(t) (short dot) and 044 (i) 
(dash-dot-dot) are presented for comparison as well. 

Figure 6. The logarithms of the relative percentage errors for the Fade approximant 
Pi{x) (dash) and for the super-exponential approximant (pl^it) (solid), in the case of 0(i) = 
sin(t)2exp(-t2). 

Figure 7. Multipliers M*2{t) (dash-dot), M;^{t) (dot), M^(i) (dash), M*^{t) (solid 
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line) in the case of = ln(l + t) exp(— i^) . 

Figure 8. The self-similar exponential approximants 044(i) (dashed line) and (j)l^{t) (solid 
line), bracket the exact function (f){t) = ln(l + t) exp(— t^), shown with dash-dot-dot. The 
approximants 022 (^) (dash-dot line) and (pl^i^) (dotted line) are presented as well. 

Figure 9. The Fade approximants Fl{t) (dash), Pl{t) (sohd), Fj{t) (dash-dot), P|(i) 
(dot) and 0|5(t) (short dot) are compared to each other and with the exact (j){t) = ln(l -|- 
t)exp(-t^) (dash-dot-dot). 

Figure 10. Logarithm of the relative percentage errors for the self-similar exponential 
approximant 055(i) (sohd line) and for the Fade approximant Pf{t) (dash), in the case of 
(l){t) = ln(l + t)exp(-t2). 

Figure 11. The dependence of the multipliers M^2{t) (dash-dot), Ml-^{t) (dot), M^^{t) 
(dash), M^^lt) (solid )as a function of time, in the case of 0(t) — exp(— t) 

Figure 12. The approximants (f)l^{t) (dot), 044(i) (dash), (f)l^{t) (solid) are compared to 
the exact function 0(t) = j^ex.p{—t) (dash-dot). 

Figure 13. The relative percentage errors for the Fade approximant Pi(a:) (dash) and 
for the self-similar exponential approximant (pl^it) (solid) are shown in the case of 4>{t) — 
^exp(-i). 

Figure 14. Dependence of the multipliers Ml2{t) (dash-dot), Mlj^{t) (dot), Ml^{t) (dash), 
Ml^ii) (solid line) as a function of time in the case of 0(t) = exp(— t (1 + t)"^^^). 

Figure 15. The self-similar exponential approximants 044(t) (dashed line) and 055(t) 
(solid line), bracket the exact function 0(i) = exp(— i (1 -|- i)~^^^) shown with dash-dot-dot, 
rather poorly. Approximants 022(0 (dash-dot line) and 033(t) (dotted hue) are presented as 
well. 
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Figure 16. Error of the approximant (pl^it) (solid line) in the case of = 
exp(— t (1 + compared with that of the Fade approximant -P4 (t), shown with dashed 

line. 

Figure 17. Multipliers M^^it) (dash-dot), M^^{t) (dot), Ml^{t) (dash), M^^{t) (solid 
line) in the case of the Stieltjes function. 

Figure 18. Approximants (l)l^{t) (solid hne) and 044(t) (dashed hne) to the Stieltjes func- 
tion, are compared with the exact expression (dash-dot line) and with the Fade approximant 
Pi{t), shown with dotted hne. 

Figure 19. In the case of Stieltjes function, it is shown here that the Fade approximant 
P^{t) shown with dashed line, easily outperforms (pl^it) shown with solid line 

20. Illustration of how well the "exact" numerical solution to Eq. (11), can be ap- 
proximated by the superexponential approximants in dimensionless units X and T , for 
X(0) = 50, m = 0.85. The "Exact" numerical solution to the Eq. (11) is shown as the 
solid line; the second -order approximation (15) is shown as the dash-dot-dot line; the third 
order approximation is shown as the dash-dot line; the fourth order approximation is shown 
with dotted line, while the best, fifth-order approximation, is shown with the dashed line 
and can barely be distinguished from the exact solution. The naive exponential expression 
xjv(i) (13) is shown with the short dash line. 

Figure 21. Demonstration that the two-point Fade approximant Pi{t) (16) shown with 
dashed line, completely fails at intermediate times. The approximate solution (15), based 
on the fifth-order superexponential approximant and shown with solid line, works well. 
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